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We propose a method of diagnostic of a degenerate ground state of Bose condensate in a double 
well potential. The method is based on the study of the one-particle coherent tunneling under 
switching the time-dependent weak Josephson coupling between the wells. We obtain a simple 
expression that allows to determine the phase of the condensate and the total number of the particles 
in the condensate from the relative number of the particles in two wells An = m — ni measured before 
the Josephson coupling is switched on and after it is switched off. The specifics of the application 
of the method in the cases of the external and the internal Josephson effect are discussed. 
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PACS numbers: 03.75.Fi, 05.30.Jp 

Beginning from its first observation [Ij the Bose- 
Einstein condensation (BEC) of atoms in alkali metal 
vapours remains the source for new possibilities for the 
study of macroscopic quantum phenomena. One of these 
phenomena is the coherent tunneling of atoms between 
two coupled Bose condensates (BC) @, that is analogous 
to the Josephson effect in superconductors. It is known 
[3j that for the case when the total number of atoms in 
the trap TV = n\ + n 2 is conserved and the trap is sym- 
metric relative two BC the average relative number of 
atoms ri\— ni = (^Ini — J^l 1 ^) is equal to zero in the 
ground state and in any excited state. Therefore, one 
can expect that in such a situation the study of nonsta- 
tionary coherent tunneling (which is realized when one 
or several parameters of the system depends on time) is 
more informative for the diagnostics of the macroscopic 
wave functions of the condensates. In this case the av- 
erage value of the relative number of the atoms n\ — n 2 
measured in the time to is generally nonzero one and it 
depends on the history of the systems at all t < t . In 
this paper we show that nonstationary Josephson effect 
can be used for the diagnostics of a macroscopic state of 
BC and the total number of the atoms in the condensate. 

We consider the simple model of coherent tunneling 
between two BC, described in Q| (see also references in 
it). The model is based on the two mode approxima- 
tion, that implies that each of N bosons can be in one 
of two states and the dynamical coupling between these 
state allows the bosons to jump from one state to the 
other. Such a model is applicable for a description of the 
external as well as the internal Josephson effect in Bose 
systems. The external Josephson effect |5| can be real- 
ized if the Bose gas is confined in a double- well trap and 
the tunneling between two wells is small. In this case two 
modes correspond to self consistent ground states in each 
well. The internal Josephson effect [6| can be realized in 
a Bose gas with two macroscopically occupied hypcrnnc 



states (e.g. the \F — I, nip — —1) and \F — 2, nip = 1) 
states of 87 Rb atoms). The dynamical coupling between 
the two states is settled by a resonant laser field applied 
to the system. At the beginning we specify the simplest 
case of the external Josephson effect at T = 0. 

The Hamiltonian of the symmetric two-mode model 
has the form 



H = H°+H 1 {t) 



^(m-n 2 r-^-(ata 2 +h.c), (1) 
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where aj~ (a,) are the creation (annihilation) operators 
for the well i, and hi = a^ai are the number operators. 
The parameter K is determined by the interaction be- 
tween the atoms in the well. Here we consider the case 
of K > that corresponds to the repulsive interaction. 
The value of Ej is determined by the overlapping of the 
wave functions of two modes and it can be controlled by a 
variation of the height and (or) width of the barrier. For 
the external Josephson effect without loss of generality 
one can set the Josephson coupling Ej(t) be real. 

Let us consider the situation when the dynamical cou- 
pling between two condensates is switched on at t = 
ti = and switched off at t = tf. At t < and 
t > tf the coupling parameter Ej(t) = and the oc- 
cupation numbers operators h\ and h 2 as well as the rel- 
ative number operator h\ — h 2 commute with the Hamil- 
tonian and do not depend on time. During the time 
when the coupling is switched on the operator h\ — hi 
is changed. Let at t = the wave function of the two 
mode Bose condensate is ^(0) and at t = tf it be- 
comes ^{tf). The task we consider is how to find the 
characteristics of the function ^(O) from the measure- 
ments of the mean relative number (^')|ni — hi\ty). Let 
us specify the case of an odd total number of the par- 
ticles (the case of an even N is discussed below). At 
N = 2A1 + 1 and Ej = the ground state of the Hamil- 
tonian J^ is double-degenerate. The minimum of the 
energy equal to K/8 is reached for the orthogonal states 
\gi) = \M + 1,M) and \g 2 ) = \M, M + l) as well as for an 
arbitrary superposition of these states \g) = a\gi) +b\gi) 
(\a\ 2 + \b\ 2 = 1). The state at t = can be parameterized 



as |tf(0)) = cos(0/2)|M + l,M) + sm(0/2)e^|M,M+l) 
At 9 ^ 0, 7r this is the entangled state. The angle 9 is 
connected with the initial relative number by the relation 



An(0) = <tf(0)|ni - fta|*(0)> = 



cos( 



(2) 



Since this value does not depend on tp the phase can- 
not be extracted from the result of measurements of the 
initial relative number. But the phase tp is also an es- 
sential characteristics of the macroscopic state of the 
Bose condensate. In particular, the interference pat- 
tern emerging under an overlapping of two BC with the 
phases tpi and tp 2 is determined by its relative phase 
4> = if i — <f2- We will show that the value of tp can 
be determined from the measurements of the final rela- 
tive number An/ = (^f(tf)\ht — fi2\^{tf)). To do this 
the amplitude of the Josephson coupling should be taken 
such a small that the strong inequality NE 1 J lax / K <C 1 
be satisfied. Then at < t < tf the system remains 
in the Fock regime. In this regime the dynamics of 
the system is realized mainly on the states for which 
|(*(*)|rii - h 2 \ *(*))| < 1. Therefore to find the evo- 
lution of the function *f> one can use the basis (\gi), |<?2))- 
Note that the reqime considered is the same as required 
for a realization of the Bose qubit Q • 

It is more convenient to use the unitary transformed 
basis of symmetric \s) = (\gi) + | (72) )/v / 2 and antisym- 
metric I a) = i\g\) — \g%)) / \2 states. In this basis the wave 
function of the BC reads as W(i) = s(t)|s) +a(t)\a) Using 
the nonstationary Schrodinger equation ifi& — H 1 ^ one 
finds that the functions s(t) and a(t) satisfy the equations 



ihs 



K Ej{t)(M- 



1 



-s, 



n- K 

ina = — a - 



E 3 (t){M + l) 



(3) 



(4) 



Integrating equations @, j|| we find the mean value of 
the relative number at the time tf 

An{t f ) = s*(tf)a(t f ) + c.c. = s* (0)a(0)E- 2iQ + c.c. = 

= cos 9 cos(2$7) — sin 9 sin(2Q) sin tp,(5) 

where ft = (l/2K)(M + 1) f s Ej(t')dt'. 

Eq. (JSJ determines the relation between the measured 
quantity An(£/) and the parameter tp. One can see that 
for entangled initial states the relative number in a final 
state depends on the phase tp and this phase can be found 
from the measurement of An(£/). 

Thus, if one has a system in a reproducible (but un- 
known) initial state |^(0)) the parameters 9 and tp that 
describe this state can be found from two sets of mea- 
surements of the relative number at t = and t > tf 
(under assumption that the total number of particles N 
in the condensate is known) . If the total number of the 
particles is unknown an additional set of measurements 
is required: the measurement of the final relative number 
occupation with another value of tf . Using the results of 



three sets of measurements one can determine the initial 
state and find the total number of the particles in the 
condensate. 

It is necessary to point out an essential restriction for 
the maximum value of tf — U. In deriving JSJ) we do 
not take into account that the coupling between two 
condensates induces small (of order of Ej(M + 1)/K) 
but nonzero occupation of the excited states \e\) = 
\M + k + l, M - k) |e|) = \M - k, M + k + 1) (with 
k > 0). Due to such processes the phases of s(t) and a(t) 
are shifted from the values given by the solution of Eqs. 
© , 0) ■ If such a shift is of order of unity the relation (JSJ 
is not valid any more. Nevertheless, one can show that 
for tf-ti < hK/[EJ ax {M+l)f the phase shifts are very 
small and Eq. JSJ is applicable. The fulfillment of the 
mentioned restriction on the value of tf — ti is required 
for the use of the method of the diagnostics proposed. 

Let us now discuss the case of the BC with an even 
number of atoms. In the symmetric double-well trap 
the ground state of the condensate with N = 2M 
is \g) — \M,M). This state is nondegenerate and 
An(0) = An(tf) = 0. If initially the system in an ex- 
cited state |e) = a\M + 1, M - 1) + b\M - 1, M+ 1) then 
An(0) = 2(|a| 2 — \b\ 2 ) can be nonzero, but for tf — U <C 
hK/[EJ ax (M + l)] 2 the difference An(t f ) - An(0) is of 
order of ME™ ax j K ^C 1 Such behavior are in difference 
with the case of odd N when the change of An can be 
order of unity. This feature can be used for determining 
the parity of the number of atoms in BC. We point out 
again that this conclusion is for the symmetric relative 
two BC confining potential. 

If the confining potential is an asymmetric one the 
Hamiltonian Q is modified to 

— -s-("i ~ ^a) 2 - A/j,(hi - n 2 ) — (afa 2 + h.c).(6) 

o Z 

One can see that if A/i = K/A : the ground state of the 
system with an even number of the atoms N = 2M is 
double-degenerate at Ej = 0, and its wave function can 
be presented in the form \g) = a\M+l, M-l)+b\M, M). 
This situation is analogous to the symmetric case with 
odd N. The only difference is that the values of An(0) 
and An(tf) given above are counted from An = 1. Thus, 
under assumption that one can control the value of A/i 
with the accuracy |A/i - K/4\ < MEf ax the method 
of diagnostic of the ground state wave function and the 
total number of the atoms suggested is applicable for the 
BC with even N. 

Hitherto we have discussed the case of the external 
Josephson effect. The case of the internal Josephson ef- 
fect is also described by the equation © (in the rotating 
frame of reference) |j] . In this case the expression for the 
chemical potential A/i reads as 

. ANirh 2 . 

Afj, = -d-\ 77(aii-a 22 ), (7) 

m 

where S is the detuning of the laser field from the resonant 
frequency, an and 022, the s-wave scattering amplitudes 



of macroscopically occupied internal states |1) and |2), 
correspondingly. In a situation, where the value of the 
detuning can be varied smoothly, one can achieve the 
regime of the degenerate ground state of the Hamiltonian 
Hq both for the even (6 = 6 e ) and for the odd (S = 5 ) 
number of atoms N. In such a regime one can apply the 
method of diagnostics of the initial state of BC proposed 
in this paper. We note that for the case of the internal 
Josephson effect the value of An(tf) is just proportional 



to the expectation value of M z (tf) - the projection of 
the magnetic momentum of the BC on the axis chosen. 
Therefore, using Eq.© one can determine the phase <p 
and the total number of the atoms in the condensate from 
the measurement of M z (tf ) . 
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